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Abstract 

We report a measurement of the forward-backward charge asymmetry of electrons from W boson 
decays in pp collisions at yfs = 1.96 TeV using a data sample of 170 pb _1 collected by the Collider 
Detector at Fermilab. The asymmetry is measured as a function of electron rapidity and transverse 
energy and provides new input on the momentum fraction dependence of the u and d quark parton 
distribution functions within the proton. 

PACS numbers: 13.38.Be, 13.85. Qk, 14.60. Cd, 14.70.Fm, 
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I. INTRODUCTION 



A necessary input for cross section calculations at a hadron collider is an estimate of 
the momentum distribution of the incoming partons that participate in the hard-scattering 
process. The probability of finding a parton carrying momentum fraction x within the 
incoming hadron is expressed in the parton distribution function (PDF). At the Tevatron, 
any cross section calculation will have to integrate over the proton and anti-proton PDFs. 
Presently, many measurements at the Tevatron have significant uncertainties associated with 
the choice of PDF. These uncertainties will become more important as the datasets continue 
to grow. For example, PDF uncertainty is expected to be among the dominant systematic 
uncertainties in a precision determination of the W boson mass. 

The PDFs are not calculable and must be determined using measurements from a wide 
range of scattering; processes . Measurement of the forward-backward charge asymmetry 
in pp —>■ + X provides important input on the ratio of the u and d quark components 
of the PDF. Since u quarks carry, on average, a higher fraction of the proton momentum 
than d quarks 3], a W + produced by ud — > W + tends to be boosted forward, in the 
proton direction. Similarly, a W~ tends to be boosted backward. This results in a non-zero 
forward-backward charge asymmetry defined as 

_ da(W + )/dy w - da(W-)/dy w 
[VW) da(W+)/dy w + da(W~)/dy w ' 1 j 

where yw is the rapidity of the W bosons and da(W ± ) / dyw is the differential cross section 
for W + or W~ boson production. 

Leptonic decays of the W boson, in our case W — > ev, provide a high purity sample for 
measuring this asymmetry. However, because pz of the neutrino is unmeasured, yw is not 
directly determined, and we instead measure 

= da(e + )/dr] e -da(e-)/dr] e 
' ' da(e + )/drj e + da(e~)/dr] e ' 



where r\ e is the electron pseudorapidity |4jj. By assuming the W — ► eu decays are described 
by the Standard Model V — A couplings, Aiiy) probes the PDF. 

Previous measurements of the asymmetry f|, using 110 pb _1 of pp data at <Js =1.8 TeV 
collected by the Collider Detector at Fermilab (CDF), have provided constraints on the 
PDFs for u and d quarks at momentum transfer of Q 2 ~ M\y 2 . In this article we describe 
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a new measurement based on data collected with the CDF II detector at y 7 ? = 1.96 TeV 
corresponding to an integrated luminosity of 170 pb _1 . We measure the asymmetry in two 
regions of electron Et that probe different ranges of yw and thus increase sensitivity to the 
PDFs in the region x > 0.3 where currently they are least constrained. 

II. DETECTOR DESCRIPTION 

The CDF II detector |f| has undergone a major upgrade since the previous data-taking 
period. The components relevant to this measurement are described here. 

Tracking detectors immersed within a 1.4 T solenoidal magnetic field are used to re- 
construct the trajectories (tracks) and measure the momentum of charged particles. The 
Central Outer Tracker (COT) is a 3.1 m long open-cell drift chamber which provides track 
measurements (hits) in 96 layers in the radial range 40 cm < r < 137 cm |2[. Closer to the 
beam, a silicon tracking system .8;] provides precise hits from eight layers of sensors spanning 
1.3 cm < r < 28 cm and extending up to 1.8 m along the beam line. The COT allows track 
reconstruction in the range |r/|^l. The silicon detector extends that range to |ry|^2.5. 

Segmented electromagnetic (EM) and hadronic calorimeters surround the tracking sys- 
tem and measure the energy of particles The energy of electrons is measured by lead- 
scintillator sampling calorimeters. In the central region, \rj\ < 1.1, the calorimeters are 
arranged in a projective barrel geometry and measure EM energy with a resolution of 
[a{E T )/E T ] 2 = (13.5%)7£ T (GeV) + (2%) 2 . In the forward region, 1.2 < \r}\ < 3.5, the 
calorimeters are arranged in a projective "end-plug" geometry and measure EM energy with 
a resolution of [a(E T )/E T ] 2 = (14.4%) 2 / E T (GeV) + (0.7%) 2 . 

Both central and forward EM calorimeters are instrumented with finely segmented detec- 
tors which measure shower position at a depth where energy deposition by a typical shower 
reaches its maximum. In the central region we use proportional wire chambers with cathode 
strip readout, in the forward region shower position is measured by two layers of 5 mm wide 
scintillating strips with a stereo angle of 45 degrees between them. 
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III. DATA SETS AND SELECTION 



Our signal sample is comprised ofW^-eu candidate events, and a sample of Z° — > e + e~ 
candidate events is used to calibrate the charge identification. Events of interest are initially 
selected by an online trigger system with differing requirements for the central and forward 
regions. For W candidates, the central trigger requires an EM energy cluster with Et > 
18 GeV and a matching track with pt > 9 GeV/c. To avoid any potential charge bias in the 
track trigger efficiency, we also accept events from a trigger which requires an EM energy 
cluster with Et > 20 GeV and missing transverse energy (Et ) of at least 25 GeV, but has 
no explicit track requirement. The forward trigger for W candidates requires an EM energy 
cluster with Ex > 20 GeV and $t > 15 GeV. A backup trigger drops the $t requirement 
and is used to estimate the QCD jet background contribution. The trigger for Z candidates 
requires two EM energy clusters with Et > 18 GeV. 

The criteria used to identify the electron and positron candidates, which are described 
in detail in reference nj] and summarized below, are designed to reject the energy deposits 
from photons or QCD jets. 

• Et > 25 GeV, 

• F Iso < 0.1, where F Iso = additional energy in an "isolation" cone, of angular radius 



R = a/ (A0) 2 + (A77) 2 = 0.4 centered on the electron, divided by the electron energy, 

• A small amount of associated hadronic energy, less than 5% of the EM energy, 

• The shower shape in the EM calorimeter and shower maximum detector must be con- 
sistent with that observed from test-beam data, 

• The position along the beamline of the pp collision, z , is well reconstructed and \z \ < 
60 cm H, 

• A track consistent with the position and energy measured in the calorimeter. 

COT tracks, reconstructed independent of the calorimeter measurement, can be compared 
to it in position and momentum. However, the coverage of the COT is limited to |//| ~ 1. To 
extend the measurement to higher \rj\, we instead use silicon tracks reconstructed by a new 
calorimeter-seeded algorithm as described below. Two points and a signed curvature define 
a unique helix. The positions of the electromagnetic shower and of the pp collision provide 
the two points. The curvature of the trajectory is predicted from the transverse energy 



measured by the calorimeter. These two points and the curvature are used to generate two 
seed helices and associated covariance matrices, one for each charge hypothesis. Those seed 
helices are then projected into the silicon detector where hits are attached using a road- 
based search and requiring at least 4 attached hits with x 2 /dof < 8. If silicon tracks are fit 
for both charge hypotheses, the x 2 /dof is used to identify the charge with the best fit, and 
cases with A% 2 /dof < 0.5 are rejected as ambiguous. 

The relative alignment of the silicon detector and the calorimeter is determined using a 
sample of well identified with both COT and silicon tracks. To avoid a charge bias from 
the W charge asymmetry, we explicitly equalize the number of events of each charge used 
in the alignment for rj > and separately for rj < 0. Offsets of (9(1 mm) and rotations of 
O(10 mrad) are measured and corrected. The resulting position resolution in the forward 
calorimeter is measured to be 1 mrad, consistent with the design expectation. 

Candidate W — > ev events are required to have exactly one such e ± candidate as well 
as $t > 25 GeV and transverse mass in the range 50GeV/c 2 < Mt < 100GeV/c 2 . To 
suppress backgrounds from QCD and Drell-Yan processes, we require that there be no other 
EM energy depositions with Et > 25 GeV. The selected sample contains 49,124 central and 
28,806 forward events. 

IV. MEASUREMENT OF THE CHARGE ASYMMETRY 

Directly measured in the experiment and shown in Figure ^ is the raw, uncorrected, 
asymmetry. In order to reconstruct A(rj e ), the measurement needs to be corrected for 
the effects of charge misidentification and background contributions. These rj dependent 
corrections are applied bin-by-bin, and binning coarser than shown in Figure ^ is used to 
reduce the effect of the uncertainty from these corrections. 

A. Charge Misidentification 

The electron identification is constructed, and observed, to have a charge symmetric 
efficiency However, resolution effects can lead to misidentification of the charge, which 
dilutes the asymmetry. Residual misalignments in the silicon detector and calorimeters could 
give rise to a bias in the charge identification that would directly bias the asymmetry. We 
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FIG. 1: The raw, uncorrected, charge asymmetry is plotted as a function of electron rj. 

measure the probability of such misidentification and correct for it. Calling that probability 
/+ for e + and /_ for e~, the corrected asymmetry can be computed from the raw asymmetry 
as A = (A raw + /+ - /_)/(l -/+- /_). 

We measure f±{i]) with Z° — ► e + e~ events where a track matched to one lepton tags the 
charge of the other. The tagging leg must have \rj\ < 1.5, and COT track information is 
used if it is available. The average misidentification probability i.e., without distinguishing 
between e + and e~, is shown as a function of r\ in Figure El The difference between the 
misidentification probability for e + and e~ is shown in Figure El 



B. Background Corrections 

We correct the measurement for the contributions of three sources of background: QCD 
jets, Z° — > e + e~, and W — ^ ru — > evvv. 

The background contribution from QCD jets faking the W — > ev si gna ture is measured 



by comparing the isolation of the e ± candidate to the $t m the event 10(. Electrons from 
W decays tend to be isolated i.e., have low F iso values, while background from QCD jets 
have larger values. Similarly, W — > ev events have large ^ while QCD jets have lower 
values. If there is no correlation between isolation and $t for QCD jets, we can measure 
their shapes in the non-IV regions and extrapolate them into the signal region. Studies of 
these variables demonstrate that they are not correlated if the selection requirements related 
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FIG. 2: The charge misidentification probability is plotted as a function of electron 77. 
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FIG. 3: The difference in charge misidentification probability of e + and e , f+(rj) — f—(jj), is 
plotted CIS Si function of electron 77. 

to the EM shower shape are relaxed. Including those requirements suppresses events with 
high values of _Fj so , which makes the extrapolation statistically imprecise and degrades our 
ability to measure any potential correlation, so we remove them in estimating the QCD jet 
background. That results in an overestimate of the background, but it yields a statistically 
and systematically robust estimate. This measured upper bound on the background fraction 
is 2% for \rj\ < 1 and increases to about 15% for \rj\ > 2. We correct the raw asymmetry by 
a factor of 1 + -Fqcd, where for the background fraction, Fqcd, we use half the measured 
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upper bound, with uncertainties of ±50%. Since we have only an upper limit, this choice 
provides full coverage of the actual value at 2a. 

Z° — > e + e~ events in which one of the leptons is lost represent a small, but asymmetric 
background 1 121. This background contribution is determined with a Monte Carlo calculation 

n 

using the PYTHIA generator [131] . and it corresponds to about 1% of the signal. W — > ru — > 
evvv events bias the measured asymmetry because the r decay dilutes the information 
available in the direction. This background contribution is about 4% of the signal. 
The number of e + and e~ events predicted for these backgrounds are subtracted from the 
measured values bin-by-bin in rj. 

Figure |U shows the fully corrected A(r] e ). 
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FIG. 4: The fully corrected charge asymmetry is plotted as a function of electron rj. Both statistical 
and total (statistical+systematic) uncertainties are shown. 



C. Et dependence 

The asymmetry probes a large range of x for the parent u and d quarks, from an upper 
value of approximately 0.5, where valence quarks dominate, down to 2 x 10~ 3 , where sea 
quarks dominate. Large values of yw correspond to the extreme values of x. For example, a 
high-x u quark and a low-x d quark lead to W + with large pz and therefore large yw The 
V — A couplings in the W — > ev decay cause the e + to be preferentially emitted opposite the 
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W + flight direction. The electron asymmetry, A(j] e ), is a convolution of these competing 
production and decay asymmetries, which results in the sign change of A(i] e ) at large |?7 e |. 

Direct sensitivity to the PDF would be improved by reducing the decay asymmetry effect 
e.g., by reconstructing the W direction. The unmeasured pz of the neutrino and the poor 
$t resolution complicate this reconstruction. However, we can improve the correspondence 
between r] e and yw based on the kinematics of just the electron, which is well measured. The 
neutrino pz ambiguity is a smaller effect for electrons with high Et than for those at low Et- 
We exploit this by separating the asymmetry measurement into bins of electron Et- The 
size of the statistical and systematic uncertainties allow two bins, 25 GeV < Et < 35 GeV 
and 35 GeV < Et < 45 GeV. For a given r] e , the two Et regions probe different ranges 
of yw, and therefore x, and the higher Et bin corresponds to a narrower range. As a 
result, measuring the asymmetry separately in the two bins allows a finer probe of the x 
dependence. 

D. Systematic uncertainties 

The corrections for charge misidentification and background contributions are measured 
and applied separately for each Et bin since they are Et dependent. The statistical uncer- 
tainty on the charge misidentification correction dominates the systematic uncertainty on 
the asymmetry measurement. The uncertainty from the QCD jet background correction is 
small, and the other background uncertainties are negligible. 

Detector misalignments can induce an inherent charge bias. Such biases would be 
naturally corrected by the charge misidentification probabilities measured from the data. 
Nonetheless, we check the robustness of the charge determination by varying the alignment 
corrections within their uncertainties and verifying that the resulting changes in the asym- 
metry are not significant. We also verify that using COT tracks, when they are available, 
instead of silicon tracks results in no significant difference. 

CP invariance requires A(—i] e ) = —A(rj e ). The fully corrected data shown in Figure 0] 
show no evidence of CP asymmetry, the level of agreement is characterized by x 2 /dof = 
9.5/11. The ±7] e data are folded together to obtain a more precise measure of A(|?y e |). 

These results are most useful as input to future global PDF fits. Such fits use Monte Carlo 
generators without a full detector simulation. We have studied possible biases introduced 
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TABLE I: The measured asymmetry values are tabulated in percent with combined statistical and 
systematic uncertainties. The listed |r/ e | is the event weighted average. Asymmetric uncertainties 
listed for some values arise because of the Poisson and binomial statistics inherent in the event 
counting. 

by detector effects by comparing the asymmetry from a PYTHIA Monte Carlo generator 
to the fully simulated results and found no significant effects. 

E. Results 

The measured asymmetry, A(|?y e |), is listed in Table U and plotted in Figure for the two 
Et regions. Predictions from CTEQ lj and MRST .2] PDFs, which fit to previous CDF 
results j^J, are shown for comparison. Those predictions use a NLO RESBOS Monte Carlo 
calculation with soft gluon resummation to model the W pt distribution, to which they can 
be sensitive [14]. Since the previous measurements upon which these predictions are based 
are least constraining for \rj\ > 1 and do not separate the Et dependence, inclusion of our 
results will further constrain future fits and improve the predictions. 
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FIG. 5: The measured asymmetry, ^4(|r^ e |), is plotted and predictions from the CTEQ6.1M (solid) 
and MRST02 (dashed) PDFs are compared using a NLO RESBOS calculation. Both statistical 
and total (statistical+systematic) uncertainties are shown. The upper plot is for 25 < Et < 35 
GeV. The lower plot is for 35 < E T < 45 GeV. 
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